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ABSTRACT 

The analysis of a proposed new technique for 
capacitance type blade tip clearance measurement 
is presented. The capacitance between the blade 
tip and a mounted capacitance electrode within a 
guard ring forms one of the feedback elements of 
a high speed operational amplifier. The dif- 
ferential equation governing the operational 
amplifier circuit is formulated and solved for two 
types of inputs to the amplifier - a constant 
voltage and a ramp. The resultant solution shows 
an output that contains a term that is propor- 
tional to the derivative of the product of the 
input voltage and. the time constant of the feed- 
back network. The blade tip clearance capacitance 
is obtained by subtracting the output of a balanc- 
ing reference channel followed by integration. 

The proposed sampled data algorithm corrects for 
environmental effects and varying rotor speeds 
on-line, making the system suitable for turbine 
instrumentation. System requirements, block dia- 
grams, and a typical application are included. 

[Key, words: Tip clearance measurement; Turbo- 

machinery; Capacitance] 

INTRODUCTION 

Computational codes for predicting gas turbine 
characteristics generated through improved graph- 
ics, numerics, and physical realism require a 
thorough understanding of the aerodynamics that 
govern combustion and heat transfer processes in 
gas turbine engines. The accuracy and predictive 
capabilities of these codes must be verified by 
comparison with reliable experimental data over 
different operating conditions. High fidelity 
instrumentation techniques that can operate in 
hostile engine environments are therefore 
required. The clearance between a turbine blade 
tip and its shroud is one of the important factors 
that govern the efficiency of a turbine engine. 
Obtaining reliable data on such clearance there- 
fore is an essential part of the above research 
objective. Earlier, different techniques were 
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used by Barranger (1,2) to obtain such data. The 
present paper gives the details of the analysis 
of a proposed system that is expected to meet the 
objectives. 

SYSTEM REQUIREMENTS 

Figure 1 is a schematic of a rotor blade and the 
probe tip of a typical shroud mounted capacitance 
sensing probe. As the blade passes under the 
probe tip, the capacitance Cp(t) between the 
inner electrode and the blade tip changes as a 
function of the geometry of the gap. If C 0 
represents any static capacitance that is present 
or that is deliberately introduced across Cb(t), 
then the capacitance at any moment is given by 

C(t) = C Q + C b (t) (1) 

Figure 2 shows C(t) as a function of time for the 
passing of two blades under the probe tip. 

The probe tip diameter usually is made much 
smaller than the inner blade spacing, which means 
C(t) will be equal to C 0 between the blades. 

The blade capacitance Cp is the difference, 
between the peak value of C(t) and Cp. This 
difference increases with decreasing clearance. 

The instrumentation system should be capable of 
measuring Cp accurately at all rotor speeds. 
Furthermore, it should also have static calibra- 
tion capability, good dynamic response, and suf- 
ficient sensitivity. The proposed system is 
expected to meet these requirements and be suit- 
able for turbine instrumentation. 

PRINCIPLE OF OPERATION 

The objective of the proposed system is to measure 
the capacitance between the blade tip and an 
electrode in a transducer mounted on the shroud. 
This electrode is surrounded by a guard ring and 
an earthed sleeve as shown in Figure 3. The 
transducer is connected to an operational ampli- 
fier circuit as shown in Figure 4. Doebelin (3) 
describes a similar capacitance feedback technique 
except that now a feedback resistor R 2 is 
introduced and a different type of input voltage 
is applied. In the present system, the total 
capacitance across the feedback resistor changes 
with tip clearance, effectively introducing a 
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change in the RC time constant of the circuit. 

This change in the tjme constant .-is measured to 
obtain the capacitance change due to tip clearance 
by the methods described below. Another feature 
of the circuit in Figure 4 is that the connecting 
cable and the dielectric material used in the 
transducer appear effectively across the two 
inputs of the operational amplifier. 

In an ideal operational amplifier the potential 
difference between its inputs approaches zero*. 
Hence, any variation in the dielectric material 
and cables should have minimum effect on the 
system performance. Further, because of the 
almost zero potential difference between the two 
inputs, the main electrode and the guard ring will 
be almost at the same potential, thus providing 
the required guarding. Various conditions under 
which the change in time constant can be measured 
are analyzed below. 

ANALYSIS 

To formulate the governing differential equation, 
the blade tip .clearance capacitance C b (t) is 
lumped with the. fixed capacitance Cq to form _ 
C(t) =. C 0 .+ Ck(t) as shown in Figure b. The oper- 
ational amplifier . is assumed to be ideal. Then, 
the current through R ^(=V i/R i ) is equal and oppo- 
site in sign to sum of the currents through R2 
and* C(t). FromFigure 5 , the current through R2 
will, be /'Vi (0/^2 and the current through C(t) 
will be d(V2 • C(t))/dt. The resulting differ- 
ential equation for the above circuit is 


■ *■ V 2 (t) = eft)* V 2 (0) 

This particular case may be difficult to implement 
in practice since to apply this method the inte- 
grator must be brought to the same initial voltage 
condition between blade pairs. Furthermore, 
static calibration is not possible, because the 
amplifier will be driven to saturation by the 
amplifier offsets. 

In cases 2 and 3 , constraints are placed on the 
R2C ( t ) time constant. As seen in Figure 2 , C b 
(corresponding to the clearance) occurs somewhere 
in the middle of the total blade passage time T. 
The R 2C ( t ) time constant is chosen to be small 
compared to T in order that the transients 
resulting from a change in the input die down well 
ahead of this peak. For a good high frequency 
response, : R2C(t) is also required to be much 
less than the period of the highest significant 
frequency in the blade tip capacitance waveform. 

CASE 2 : Vf IS A dc VOLTAGE 

The approximate solution of Equation ( 2 ) when Yj. 
is a constant voltage is given by, (refer to the . 
appendix for the derivation) 

' - v 2 (t) = 7q-[ v i - -jt ( v i K 2 c(t »] 

or 


V x (t) 

~*r 


v ? (t) dem dV 2 (t) 

+ v,(t) • + C(t) ' 4 


~ 3 F 


V 1 R 2 r 

-V 2 ( t ) = ^[1 - R; 


dc b (t) l 
' — JF~J 


( 3 ) 


or 



( 2 ) 


The objective of the analysis, is to establish a 
relationship between V2U) and C b (t) for a given 
type of input'. V].. Different configurations can 
be visualized satisfying Equation ( 2 ) depending 
upon the choice of R2 and type of V]_.. However, 
only three cases are analyzed because of their 
direct application to the blade tip clearance 
capacitance measurement. 

CASE 1 : *' INTEGRATOR WITH AN INITIAL' CONDITION ", . 

With R2 removed from ..the .circuit (R2 = . 
infinity), V^ = 0 and V 2 ( 0 ) as the initial 
voltage, the circuit becomes an integrator with 
an initial condition V2 ( 0 ). . Equation ( 2 ) . will 
then be, 

dV 2 (t) V 2 (t) : bc(t )• " ’ ■ 

- “ - dt , C(t ) dt 

As shown in the appendix, the solution of. this 
equation is, • 


In order to find C b , a solution of Equation ( 3 ) 
is required. Figure 6 is an. implementation of 
that solution. For the purpose of deleting the 
constant term in Equation ( 3 ) and thereby obtain- 
ing a voltage proportional to dC b (t)/dt, a ref- 
erence channel containing R2 and C 0 identical 
with that of the signal channel is introduced;.' 1 * 
(A means -for correcting the output when the signal 
and reference channels are not exactly identical 
is shown later.) The output of the signal channel 
is subtracted from the ref erence channel . As a 
resultvof this balancing scheme the V^ R2/R1 
term in Equation ( 3 ) is removed and the small 
quantity proportional to dC b (t)/dt only is 
amplified by the differential amplifier. If G 
is the gain of the differential amplifier, then 
the output of the differential amplifier is given- 
by. 


GR V, - 
= R x R 2 


dC b (t) 

dt 


( 4 ) 


In Figure 6 the integrator with time constant 
R4C4 must be reset for each blade to obtain 
V 0 ( t) proportional to C b (t). The required 
sensitivity can be obtained by appropriate choice 
of G, R2, R]_» V]_, R4, and C4. 


2 



This configuration of the circuit has the advant- 
age that it requires no initialization for each 
blade and the circuit is always ready to respond 
to blade tip clearance capacitance changes. The 
only disadvantage with this method is that there 
is no response to Cp for static calibration. 

As shown below, when the input is a ramp voltage 
instead of a constant voltage, a system that 
responds to static changes in Cb can be 
realized. 

CASE 3: Vj. IS A RAMP 

For this case, is proportional to time t, 
where t starts between blade pairs and lasts the 
duration of the blade passing time T. The 
approximate solution of Equation (2) when Vi = 

Kt (K is the slope of the ramp in volts/second) 
is given by (refer to the appendix for the 
derivation) , 


' V 2 (t) = ^[ Kt “It- ( KtR 2 C(t))] (5) 

Substituting C(t) = C + C.(t) and dC(t)/dt = 
dC b (t)/dt, we obtain 0 

-V 2 (t) ^ R 2 C 0 " R 2 it ( tc b^^] 

(6) 

The schematic arrangement shown in Figure 6 is 
again used except that the input voltage is now a 
ramp. The output of the reference channel is 
changed to 



The differential amplifier output can now be 
written as, (again assuming balanced channels) 


-V d (t) = G 



dt 


( 8 ) 


In a practical situation it is difficult to match 
the R 2 C 's exactly, although Ro values can be 
matched accurately. For this case, becomes, 

GKR 9 r h 1 

-V*) =1^2 it + f V C oJ 

where R ? aC is. the difference in time constants 
R 9 C ot thft channels. The integrator output will 
theft be. 


GKR? tC.(t) 

V 0 W = -R“ * + Rt (10) ' 


where 


GKR^aC 

C 0 

p = 

The signal Cu(t) can be obtained from the meas- 
ured values of V 0 (t) using the algorithm devel- 
oped below. 

SIGNAL PROCESSING 


Equation (10) shows that the output voltage is 
proportional to the product of the instantaneous 
values of time and the blade tip clearance capac- 
itance as well as p, representing the difference 
in the initial time constants R 2 C 0 of the ref- 
erence and the signal channels. When the measure- 
ments are taken at different speeds of the rotor, 
it is clear that the tip clearance can be obtained 
only after applying corrections for the instanta- 
neous values of time. A means must also be found 
to correct for p. In the proposed system these 
difficulties are circumvented by generating a volt- 
age proportional to the rotor speed and by sam- 
pling V 0 (t) between the blades to solve for p. 

The schematic arrangement for implementing this 
scheme is shown in Figure 7. The angle clock is 
a fast updating frequency synthesizer which among 
other things will produce: 

1. A frequency proportional to the speed of the 
rotor. Let this frequency be f r . 

2. A fixed number of pulses (preset) per blade on 
the rotor at all speeds. This signal can be used 
to accurately estimate the angular position of the 
rotor. Let this number of pulses (preset) be N. 


The voltage input to the integrator one whose time 
constant is R 3 C 3 will then be. 


-V 1 = mf r 

where m is the frequency to voltage transfer 
function of the F/V converter in volts/hertz. The 
output of the integrator will be. 


1 f mf r 

V 1 = R3C3 J m * V dt = R3C3- * 4 - Kt 

( 11 ) 

Substituting for K in Equation (10) and solving 
for C b (t), we get, 


c b (t) 


R 1 R 4 C 4 R 3 C 3 

Gmf tR~ 
r c 


(V 0 (t) - pt) 


( 12 ) 


The period T in seconds between the blades (Fig- 
ure 2 ) can be related to the frequency f r by, 



(13) 
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where M is the number of blades on the rotor. 

As mentioned before, the angle clock generates a 
fixed number (N) of preset pulses per blade in the 
duration T. Therefore, at all speeds, 

T =-, N aT (14) 

where aT is the duration between the pulses. 
Hence, 


to be approximately 1 pF. For good high frequency 
response the measurement system is required to 
have an R 2 C 0 time constant of less than 1 v»s - 
By choosing R 2 = 1 Kn and C 0 = 10 pF, the time 
constant becomes 0.01 ys which easily meets this 
requirement. Further, the time constant is much 
less than 20 ys, the blade passage time T at the 
maximum rotor speed. This means that the tran- 
sient which begins at the initiation of the ramp 
will disappear well ahead of the measured capaci- 
tance peak. 


Substituting Equation (15) in the Equation (12), 
we get 


RiR,C,R,C,MN aT 

. C.(t) = 1 4 - 4 - 3 - — (Vft) - pt) (16) 

b ■ GmR^t 0 

If V Q (t) is sampled N times per blade (i.e., a 
rate of MNf r ) , then 

R,R>iC^R,C,MN 

? 

j = 1 to N. 


R 1 R 4 C 4 R ^ C ^ MN T T 

c b (j) = V (j) - pj AT (17) 

b ' ' GmRij L° J 


From Figure 2 it can be seen that = 0 at t 
aT, (C b (l) = 0, j = 1) for. large N. Equa- 
tion (17) will then be, 

. V 0 (l) = p AT 


or. 


,yn 

aT 


(18) 


Substituting for p in Equation (17), we obtain. 


Values for other circuit components were chosen 
for the measurement system. With a ramp slope of 
K = 0.5 V/ys, and a differential amplifier gain of 
100, the sensitivity becomes 100 mV for a 1 pF 
blade tip capacitance. 

The dynamic response of -the system is limited by 
the slew rate of the operational amplifiers chosen 
for the circuit. Currently available high speed 
amplifiers have slew rates greater than 100 V/us, 
which is more than adequate for the application. 
One of the most severe speed requirements is the 
resetting of the integrators between rotor blade 
pairs. A circuit of a 10 ys period ramp generator 
integrator was tested under simulated conditions 
using a high speed operational amplifier. The 
oscilloscope photograph of the test signals (Fig- 
ure 8), shows that the resetting time is less than 
1 ys, which is less than 5 percent of the blade 
passage time T (20 ys) . 

An assessment was made of the validity of the 
approximate solution for V?(t) as applied to 
the measurement system. The exact solution of 
Equation (2) was compared with the approximate 
solution of Equation (6). Equation (2) was solved 
numerically for the ramp input using Runge-Kutta 
fourth order method. The results tally with the 
solution given by Equation (6), thereby verifying 
the approximation for this typical appl ication . 

CONCLUSION 


C b (j) = s 

j = 1 to N where. 


yj) - jyi) 


RlR.CJUC.MN 

S = — -- — = a constant 

GmRo 


(19) 


( 20 ) 


Equation (19) is a normalized value for all blade 
passing frequencies, thereby eliminating the 
requirement for measurement of time t. Thus any 
change in C 0 can be corrected continuously as 
long as C 0 remains relatively constant during 
T. This will correct for drifts due to tempera- 
ture changes, other environmental effects, and 
long term aging. 


A proposed new technique for capacitance type 
blade tip clearance measurement has been 
described. An analysis was performed on an opera- 
tional amplifier with a parallel combination of 
blade tip capacitance and a resistance in the 
feedback. When the operational amplifier is 
driven with a ramp input voltage, the tip clear- 
ance capacitance can be measured with an algorithm 
which can be easily implemented. The measuring 
system has static calibration capability, is 
independent of changes in rotor speed and is not 
affected by environmentally induced changes in 
probe and circuit parameters. In a typical appli- 
cation the system has a sensitivity of 100 mV/pF 
and a time constant of 0.01 ys. 

APPENDIX 


TYPICAL APPLICATION 

Sample calculations, integrator reset tests, and 
computer simulations were performed for a typical 
application. The candidate engine has 50 blades 
on the rotor and a maximum speed of 60,000 rpm. 
The maximum blade capacitance Ct> is expected 


Assuming that the operational amplifiers used are 
ideal, the governing differential equation that 
relates the output voltage V 2 of the opera- 
tional amplifier to changes in the feedback capac- 
itance as given by Equation (2) is repeated below. 
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dV 2 (t) 

dt 


1 

= TTt7 


Vj(t) 


+ V 2^ + 


(Al) 


The circuit configuration for the above equation 
is shown in Figure 5. Let 


P = 


1 


1 . jC It) 


and 


R 2 C(ty + CTtT * dt 


V^t) 

~ Q = mtj 


(A2) 


(A3) 


The solution of Equation (Al) is known to be, 
V 2 (t) = const • e~ /Pdt + e" /Pdt J Q • e ;Pdt • dt 


(A4) 


where. 


e' Pdt = e 1 


|(R 2 TitT + -cfty ' • d 4r i ) dt 


= e 


J FrjriTy dt J" t^CTTT ’ dt 

(A5) 


-/Pdt 1 . a J "2 

e = cTtT e 


= C(t) • e 

- Joity * dt 


(A6) 


CASE 1 


With R 2 removed from the feedback (R? = infin- 
ity), V 1 = 0, and V 2 ( 0 ) as the initial voltage, 
Equation (A4) becomes. 


For the other two cases, let the time constant 
R 2 C(t) be small when compared to other periods 
of interest. In Equation (A4), this makes the 
initial transient term disappear quickly when 
compared to other term. Then, 


V 2 (t) = e ;Pdt | Q • e /Pdt • dt (A7) 


Substituting for Q, e /Pdt , e~ /Pdt , from the previ- 
ous equations, we get. 


- J R^CTtT dt r 
" V 2 (t) = “ R x C(t) * J 




city 


dt 


Rewriting Equation (A8) and integrating by parts, 


'I«^Tty dt r ( 

'J 


= R^Tity 


dt 

'fixity 


dt 

' R 2 CCt) 


du 



C dt 

3 RTCTty 


R 2 C(t)V 


RTCTty 

e R 2 |C(t) 


dV i (t) dcm 


,»-i 1 

- z ] 


where, 


R 2 e 


f dt 

J TOty 


TTtT 


•1 


C dt 

Jcnry 


c(t) • -j t - 


(A9/A10) 


dV L (t) 


+ Vi(t) • • dt (All) 

CASE 2: Vi IS A CONSTANT 

Let Vi = V. Substituting dC(t)/dt = dC^ ( t ) /dt 
in Equation (All), we get 


dt 


Z = 


f dt 

" WI r c 

R 2 Ve 6 j J‘R 2 CTtT dC b (t) 

^ J “ 


R 2 V e 




f dt 

- j OTty 


^t 


dt 


t dr 


1 


r dt 

Jjnnry 

, 2 . dt 

5 fr 2 TTty 


du 

dt 


- J R 2 C(ty 

dC.(t) R,Ve 

* R 2 C ^ * dt = cTtJ 


f dt 

J^city 



• R 2 C(t) 

dC b (t) 

dt 

dt 

d 2 C (t) 

\ 


+ b 

C(t) dt 

(A8) 

dt* 

/ 


i* 


r dt 

jR^crty 


( A12) 
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It can be shown. that the term containing the 
integral in Equation (A12) can be neglected when 
the R 2 C ( t ) is small compared to the period of 
the highest significant frequency contained in 
Cb(t). Therefore, Equation (A12) becomes, 

z • I* -at- (M3 > 


Repeating the same procedure used for constant 
input voltage and neglecting higher order terms, 
we obtain, 


Z = R^K 


C q ♦ c b (t) - 


3 W 


dc b ( 

dt 


t)“ 


(A17) 


Substituting for Z in Equation (A10) , 


-v 2 (t) = Tq- 



dC b (t) “ 
dt 


CASE 3: V x IS A RAMP 


(A14) 


The peak of the blade tip capacitance occurs 
somewhere in the middle of the total duration for 
each blade* Therefore, it is stipulated for our 
measurements, that 


2 >:> 3C o R 2 


Then, 


Substituting V* = Kt in Equation (A10), 


where, 


R 2 K j « 2 

1 = C[tj ' e 


-V 2 (t) 


( dt 
"J R,C(t)' 


Rj [ R 2 Kt " Z ] 


•1 


C dt 

jR^cTty 


dC h (t) 


+ t 


dt 


(A15) 


C(t) 


dt (A16) 


Rewriting Z and again integrating by parts, 

It dt 

nrr 


lR 2 C.(t) 


RpKe *• f 

l = r rr , • 


f dt 

jR 2 T(ty dt 
e * RxrtT 


du 


dC b (t) 


• C(t) + t — RgC(t) 


2 [ dC (t)1 

Z = R 2 K [ C o + C b (t) + t *-VJ < A18 > 

and, 

R K f dC.(tn 

-V 2 (t) = ^|t- R 2 C o - R 2 C b (t) - R 2 t -°—J 

or, 

-V t )=^-[t-R 2 C 0 -R 2 | F (tC b {t))J 

(A19) 
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Figure 1. - Schematic of rotor blade and probe 
tip of typical shroud mounted capacitance 
sensing probe. 



Figure 2. - Capacitance as a function of time for passing of 
two blades under probe tip. 
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Figure 3. - Cross-sectional view of the capacitance 
transducer. 



Figure 4. - Arrangement of the transducer in the feedback 
of an operational amplifier circuit. 
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Figure 8. - Resetting of the integrator. 
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